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Abstract 
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machinery of quantum field theory to Brownian processes. 
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1 Introduction 



The study of Brownian motion, and Markov processes based on it, has been 
a rapidly growing subject in both mathematics and physics since Einstein's 
pioneering paper [1], hnking the macroscopic heat flow to microscopic Brown- 
ian motion. The rigorous probabihty theoretical formulation is due to Wiener 
[2], leading to a well-defined Feynman-Kac path integral representation for 
solutions of the heat equation. Remarkable non-trivial results, as the non- 
differentiability of Brownian paths [3] , follow from Wiener's construction. To- 
day, Brownian motion is a very active area of research in both pure mathe- 
matics and physics [4] [5] [6], as well as in applications to a wide variety of 
phenomena in nature and finance [7]. 

Multi-particle Brownian systems with dynamical creation and annihilation of 
particles, such as branching and dying processes, are conventionally studied 
by non-linear extensions of the heat equation. Such a treatment, however, re- 
quires tailor-made kinematics for each individual system at hand, in order to 
encode an intrinsically multi-particulate problem in the formalism of a single 
heat equation. This method becomes almost prohibitively difficult in appli- 
cation to the currently intensively studied catalytic processes [8]. The aim of 
this paper is to remedy this situation and to devise a generical kinematical 
framework for any Brownian process featuring a dynamically changing parti- 
cle number. 



Our investigation starts with the crucial observation that the phase space of 
the heat equation carries a product structure, in contrast to the complex struc- 
ture of the phase space for the Schrodinger equation. The same way in which 
the latter gives rise to the ubiquitous occurrence of the complex numbers 
in quantum theory, the product structure can be absorbed into the commu- 
tative unit ring P of pseudo-complex numbers, governing the kinematics of 
Brownian particles. The study of one-particle Brownian motion then becomes 
a Hilbert module theory over P, where the non-differentiability of Brownian 
paths emerges as a direct consequence of the pseudo-complex structure of heat 
phase space. This thorough geometrical understanding of the one-particle case 
can then be put to use in the construction of multi-particle systems, employ- 
ing a Fock space approach. In contrast to the conventional treatment, the 
presented pseudo-complex Hilbert module description of the single particle 
is extended to a model-independent Fock space formalism, and the dynam- 
ics of an operator-valued second quantized field are derived. In this set-up, 
we discuss general interacting pseudo-complex quantum field theories as the 
appropriate kinematical framework for the discussion of Brownian processes 
with non-constant particle number. 
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In general, the advantage of this field theoretical point of view is that one can 
easily write down models for arbitrarily complex Brownian processes. Ideally, 
these models have exactly solvable Dyson- Schwinger equations. Otherwise, a 
truncated series expansion of interesting quantities in terms of Feynman dia- 
grams may be just the appropriate tool to extract information on analytically 
inaccessible systems. 

The abstract formalism is applied to the classically well-studied binary branch- 
ing Brownian motion. The Dyson-Schwinger equations of this system can be 
exactly solved, and coincide with the equations of motion that are found clas- 
sically for corresponding quantities. A point in case is the reproduction of the 
extinction probability for such a process, which appears as the dressed one- 
point function in the pseudo-complex quantum field theory. The formalism is 
versatile enough to embrace models with dynamical catalysts, which currently 
receive much attention [8]. 

The organization of the paper is as follows. Section II reviews the conven- 
tional treatment of multi-particle Brownian processes for the benefit of the 
non-specialist reader. Section III identifies the product structure of the heat 
phase space, and introduces the pseudo-complex ring. Section IV discusses 
a single Brownian particle in a pseudo-complex Hilbert module formalism, 
preparing the non-interacting multi-particle theory developed in Section V. 
Section VI studies theories with dynamical particle creation and annihilation. 
These abstract developments are applied to the binary branching model and 
a self-intoxicating particle in Section VII. We summarize and conclude in Sec- 
tion VIII. 



2 A review of Brownian processes 

The Brownian motion of a single particle is well-understood and can be given 
a rigorous probability theoretical description. More complex systems, such as 
spatial branching processes, are of current research interest [4] [5] [6] [9]. They 
are conventionally studied as extensions of the single-particle problem, but the 
kinematical framework has to be adapted from case to case. This section briefiy 
reviews some of the well-studied processes derived from Brownian motion for 
the non-specialist reader, and closes with a critical discussion. In particular, we 
will argue that a general framework for the study of systems with non-constant 
particle number is desirable, such that a system is completely specified by 
its dynamical equations, without the need to adapt the underlying concepts. 
The construction of such a system-independent kinematical framework for 
Brownian processes is the aim of this paper. 
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2. 1 Heat equation and free Brownian motion 



Recall that the heat equation is given by 



dt(p{t,x) = -Aip{t,x), 



(1) 



where (p is a real- valued field on non-relativistic space time ]R° x K^, and the 
diffusion constant is set to unity. Given the initial condition (/^(O.x) = m(x), 
the linear partial differential equation (1) has a unique solution. Its physical 
significance is the following: suppose the function (p{t, x) represents the tem- 
perature at time t and position x, then starting with the initial temperature 
distribution u, the heat equation describes the heat fiow in time. In this sense, 
the heat equation governs the macroscopic behaviour of the temperature field. 

In 1905, Einstein discovered the relation between the heat equation and a 
stochastic process [1] called Brownian motion. Intuitively, this process de- 
scribes the microscopic behaviour of a single particle that contributes to the 
heat. Mathematically speaking, a Brownian motion with start in a; e M*^ is a 
stochastic process, i.e., a random path {Bt : t ^ 0} which is continuous and 
has independent, stationary and normally distributed increments. A Brow- 
nian motion can therefore be considered, in a rigorous probabilistic manner 
[2], as a random element of the canonical Wiener space (C([0, oo)), JF, P^,), 
where P^; is a measure on the space (C([0, oo), JF) of continuous paths in 
equipped with the cr-field which is generated by the family of projections 
TTt : C[0, oo) — > M^, u; i— > LOt- The Wiener measure satisfies, for all times 
ti, t2,...,tn and all Borel sets Ai,A2, . . . , ^„ C R'^, the relation 



where the transition density of Brownian motion is given by the heat kernel 



¥,{Bt, e A^, Bt, eA2,...,Bt^e 






(3) 



The initial value problem for the classical heat equation. 





(4) 



is solved using the Brownian semigroup operators Pt, 
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i.e., as a linear superposition of the fundamental solutions pt- 

The semigroup property of the operator Pt, 

PtoP, = P,+„ (6) 

implies the so-called Markov property of Brownian motion, i.e., any point on 
the Brownian path can be considered as the starting point of a new Brownian 
motion. In terms of the transition density pt, the Markov property hence takes 
the form 

Pt+s{^, y) = / d'^^Pti^^ z)Ps{z, y). (7) 

By definition of the expectation, we can rewrite (5) as an integral over the 
space of continuous functions with respect to the Wiener measure P-^, 

ip{t,^)=E,{u{Bt)}= I f,{duj)u{ujt). (8) 

JCo[0,oo) 

Hence, the solution of the initial value problem for the classical heat equation 
can be represented either using the Brownian semigroup or, equivalently, as 
a path integral. An important characteristic of Brownian paths was found by 
Paley, Wiener and Zygmund [3] , who showed that with respect to the Wiener 
measure, almost all Brownian paths are nowhere differentiable. We will come 
back to this point in section 4. 

We emphasize that the heat equation (1) describes a single Brownian particle, 
by the heat function (/7(t,x). As in the case of the Schrodinger equation, the 
formalism, as it stands, cannot be extended to accommodate a non-constant 
particle number, as this would require a non-constant number of functions (fi 
in the dynamical equations. However, the desire to study killing and branching 
processes has led to the following two constructions. 



2.2 Softly killed Brownian motion 



It is interesting to study the heat equation with a potential term i'(x), which 
is assumed to be bounded from below, 

dt^{t,^) = iA(^(t,x) -^;(x)<^(t,x). (9) 

Given the initial condition (^(0,x) = m(x), the solution to (9) is given by the 
Feynman-Kac formula [7] 

(^(i,x) =E,|ii(St)exp(-^*o?s^;(5,))|. (10) 
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For a constant potential f = 7, the solutions (p{t, x) can be interpreted as 
the transition densities for a softly killed Brownian particle, surviving a time 
interval At only with probability exp(— 7At). The terminology refers to the 
fact that we are dealing with an inherent one-particle formalism, and killing 
is a matter of interpretation of the underlying random variables. This is seen 
most clearly when considering the continuous paths in Wiener space, which do 
not stop at a finite parameter value ('hard killing'), but must, in this context, 
rather be understood as random elements in C([0, 00)) x M+, i.e., as paths 
with a killing time assigned to them. 

This extension of the Wiener space as the underlying probability space is 

necessary in order to encode a process with non-constant particle number, 
in a setting that by construction can only deal with a constant number of 
particles. 



2.3 Branching Brownian motion 

The standard approach to branching Brownian motion [9] is the study of trees 
composed of independent, softly killed one-particle Brownian motions, where a 
dying particle gives birth to k children, according to a probability distribution 
{pk) on {0, 1, 2, ... }. This is formalized by consideration of a measure- valued 
process {Xt : t ^ 0}, defined by the counting measure 



where C is a Borel set in W^. Hence, if we represent each individual particle 
at a point x G M*^ by a unit point mass 5x and write for the position of the 
i-th particle alive at time t, we have 



As the exponential distribution has no memory, the measure-valued process 
X inherits the Markov property from the underlying Brownian motion. Due 
to Mc Kean (1975), the distribution of branching Brownian motion can be 
characterized in terms of a non-linear extension of the heat equation: Let ^jJ 
be a positive, bounded and continuous function on R'^ such that ^ 0(x) ^ 1 
for all X &W^. Then the quantity 



Xt{C) :— ^{particles in C at time t}, 



(11) 





(13) 
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satisfies the evolution equation 



dtip^\I^ip + ^{W{ip)-ip) ^^^^ 
(^(O,x)=0(x) 

where W{ip) := J2'kLoPk^'' is the probabihty generating function of the distri- 
bution {pk), controUing the fertihty of a dying particle. Now u is interpreted as 
describing an average mass flow of the measure- valued process X. This looks 
like a straightforward extension of the original heat equation by a non-linear 
term. However, the paths are now measure-valued, and hence we again deal 
with a somewhat different mathematical object, despite the formal similarity 
of equations (1) and (14). These adaptations were necessary to construct a 
spatially branching model within a formalism that originally described a one- 
particle system. 



One can extract the extinction probability for such a branching and dying 
process via the following construction. Denote the total mass process of the 
branching Brownian motion X by iV^ := Xt{M.'^). It is well-known (see e.g. 
[16]) that {Nt : t ^ 0} is a continuous time Galton- Watson branching process 
whose moment generating function is given by 

oo 

^e{t) = E{^^' I TVo = 1} = ^ F{Nt ^k\No^l}, (15) 

k=0 

for ^ ^ (and 0° = 1). The function (fg solves the initial value problem 



dtMt)-^{w{Mt))-Mt)), (16) 
<^e(o) = ^. 

Compared to equation (14), this evolution equation is missing the Laplace 
operator that generates the spatial evolution of the process. Hence, equation 
(16) describes the pure branching process, without the spatial evolution. For 
the special case 6^ = we have by (15) that ipo{t) = F{Nt = | iVo = 1}. The 
solution of the initial value problem (16), with ^ = 0, hence is the extinction 
probability for the total mass process in the limit t — > oo. In terms of the 
process described by (14), this corresponds to the probability that the process 
dies out eventually. We will find equation (16, = 0) as the differential form 
of the recurrence relation for the so-called one-point function in our field- 
theoretical treatment in Section 7. 
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2-4 Discussion of the traditional approach 



The one-particle Brownian motion is mathematically well-understood, and 
allows for a rigorous probability theoretical definition in terms of Wiener 
spaces. The desire to study more complex Markov processes, at the same 
level of mathematical rigor, led to the introduction of an exponential clock 
and the consideration of measure-valued paths. These generalizations of the 
one-particle case, however, come at the cost of an extension of the underlying 
probability spaces, generalizing paths C([0, oo)) to C([0, oo)) x M+, or even 
measure-valued paths, as explained in the preceding subsections. The neces- 
sity to perform these rather complicated and technical adaptations roots in 
the fact that these more complex Markov processes all feature a non-constant 
particle number, which must be encoded in a one-particle formalism, in a man- 
ner that varies from case to case. 

The problem of dealing with systems of non-constant particle number is not 
a new one, but appears naturally in relativistic quantum physics. There, the 
relativistic mass-energy equivalence E = mc^ allows for particle creation and 
annihilation, and hence is a feature of any quantum particle process at suffi- 
ciently high energy. It is clear that the Schrodinger equation. 



describing the quantum behaviour of a one-particle quantum system in a po- 
tential V , cannot feature such a dynamical creation and annihilation of parti- 
cles, due to the fixed number of functions ■0. An adequate solution was found 
in form of quantum field theories [10]. It was recognized [11] that quantum 
field theory can be formally applied to the discussion of Brownian motion, and 
vice versa, if one replaces the time t by imaginary time it. This ad-hoc adap- 
tation, however, bypasses the whole construction of a canonically quantized 
theory starting from its classical equations. There are no insights gained as 
to why one has to consider imaginary time, what the underlying geometrical 
structures are, or why there is no interference of heat in such a treatment. 
In fact, there exists no well-defined measure [12] for a Feynman path inte- 
gral which formally solves the Schrodinger equation (17). Therefore, it seems 
that using quantum field theoretical methods to describe Brownian processes 
might imply a loss in mathematical rigour. In this paper, however, we show 
that careful consideration of the structure of the heat phase space leads to a 
canonical quantization of heat, whose path integral is well-defined. 
In order to understand the occurrence of the imaginary unit in the quantum 
mechanical 'path integral', it is instructive to re- write the Lagrangian density 
for the free {v — 0) Schrodinger equation (17), 




(17) 



(18) 
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in terms of the real and imaginary parts a and (3 of the complex- valued prob- 
ability amplitude ip — a + i^, 



(19) 



The imaginary part /3 is recognized as the canonical momentum associated 
with the real part a, 

= 7^ = /3, (20) 
and the time evolution of the field phase space point (a, /9) thus reads 



d 



a 




-1 




a 






1 


2 





(21) 



revealing the complex structure of the field phase space of the Schrodinger 
equation. This complex structure is the actual reason for the imaginary unit 
entering the exponential in the quantum mechanical analogue of the Feyman- 
Kac path integral (10). It is therefore worthwhile to investigate the phase 
space structure of the heat equation in the following section. 



3 Heat phase space 



This chapter studies the geometry of the heat phase space. The central result 
is the absorption of the occurring structure into the ring of pseudo-complex 
numbers. For the discussion of Brownian motion, these play as important a 
role as the complex numbers do in quantum mechanics. 



3.1 The structure of heat phase space 



In order to identify the canonical momentum associated with the heat 0, we 
need to devise a Lagrangian density for the heat equation. This, however, 
can only be done by introduction of a second dynamical field r). Without 
such, any attempt to generate the first order time derivative in (1) from a 
Lagrangian density inevitably results in a total time derivative, and hence does 
not contribute to the Euler-Lagrange equations at all. We therefore consider 
the Lagrangian density 

£ = {dtiprj - ipdtTj) + WipWri, (22) 
yielding the Euler-Lagrange equations 
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dtV^ 



1 ^ 
1 ^ 



(23) 
(24) 



We observe that the additional dynamical field rj, required by a Lagrangian 
treatment of heat, is the canonical momentum density associated with the 
heat field (/?, as 

odt^ 

The phase space picture of the heat evolution, 



d 
di 



1 
-1 



-A 



(26) 



reveals that the geometry of the heat phase space is governed by a product 
structure [13], i.e. the representing matrix squares to the identity. The above 
matrix representation of the product structure and the phase space element 
= [(p,r)] clearly depend on a particular choice of basis in phase space. In the 
following, we will show that there is a basis-independent formulation 



Idt, 



^A0 



(27) 



for equation (26). The obvious similarity of (27) with the Schrodinger equa- 
tion (17) suggests to term the heat amplitude, and equation (27) the heat 
amphtude equation. 

Consider a linear operator / on M^, with the properties P = idR2 and / ^ 
±id]R2. Then a± = ^(idib/) are orthogonal projection operators, i.e., cr^ = a± 
and (7+(T_ = 0. As {(T+, (7_} are linearly independent over M, they provide a 
basis for R^. Thus the operator identity 



Ia± = ±(j-|- 



f28) 



can be read as an eigenvalue equation for the operator /, regarding (t± as 
vectors of M?. This shows that the matrix representation in (26) corresponds 
to the choice of the (7±-basis for (27). In particular, the heat amplitude is given 
by the R^-valued function 



ipa+ + r]a. 



(29) 



Prom this point of view, the abstract reasoning started after (26) simply re- 
produced the previous results. Now we note, however, that there is a canonical 
multiplication defined on the heat phase space R^, induced by the composition 
of operators a±. Hence, the heat phase space is not merely a vector space, but 
a commutative ring! We pause in our development of the physical aspects, in 
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order to explore this structure a bit further. 



3.2 The pseudo-complex ring 



Let F be the field of real numbers IR or complex numbers C. The pseudo- 
complex ring over F [14] is the set 

F={p = a + Ib\a,beF}, (30) 

equipped with addition and multiplication laws induced by those on F, with 
the additional rule that / ^ F is a product structure, i.e. = 1. It is easily 
verified that P is a commutative unit ring with zero divisors P'^ = P5|_ U P° , 
where 

Fl = . (31) 
It is often convenient to employ the zero-divisor representation of a number 

p = a + Ib= {a + b)a+ + {a-b)a_. (32) 

The zero-divisor branches P5j_ and P'^. are the only multiplicative ideals in P, 
thus they are both maximal ideals. Hence, the linear maps 



r± : P — ^ P/P^ ^ po ^ F, (33) 
r±(a + /6) = a±b, (34) 

are the only non-trivial ring homomorphisms from P to F, i.e., for all p,q 



r{p + q) = r{p) + r{q), (35) 

r(pg)=r(p)r(g). (36) 

Application of r± to the heat amplitude extracts the heat (p and antiheat 
77, respectively: 



r+(0) = (^, (37) 

r_(0) = r7. (38) 

Physically, this means that heat amplitudes do not interfere (linearity), and 
that one can compose heat amplitudes for complicated processes multiplica- 
tively from elementary ones. In other words, the extraction of heat from heat 
amplitudes is compatible with the Markov property (6) of Brownian motion. 
Finally, we define the linear map * : P — > P by its action *a±. = a^:. Clearly, 
* is an involution, as ** = 1. Acting on the heat amplitude, * effects an 
exchange of heat and antiheat. 
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We have seen, from the analysis of the phase space structure of the heat 
equation, that the heat amphtude is a pseudo-complex valued field. For the 
discussion of Brownian motion, pseudo-complex numbers therefore play the 
same role as complex numbers do in quantum mechanics. The systematic 
exploitation of this insight will be seen to be intimately linked to the intrinsi- 
cal properties of Brownian motion, like the non-differentiability of Brownian 
paths. 

3.3 Stationary solutions and unitary time evolution 

It is often sufficient to restrict the discussion of solutions (j){x) of the heat 
amphtude equation (27) to the solutions $(x) of the stationary heat amphtude 
equation 

^A$(x) = -£;$(x), E e R+, (39) 

which can be extended to solutions of the time-dependent dynamics by acting 
on them with the time evolution operators 

U{t) = exp{-IEt), (40) 

such that (f){t,x) — C/(i)$(x) solves (27). Note that the time evolution opera- 
tors form a semigroup, i.e., 

U{t)oU{s) =U{t + s), (41) 
and are unitary with respect to pseudo-complex conjugation, 

U*{t) = U-\t) = U{-t). (42) 

In section 4 we will identify the stationary solutions $ as elements of a Hilbert 
module 7i, and the time evolution operators (40) as unitary operators on Ti. 

3.4 Lagrangian for the heat amplitude equation 

Equipped with the constructions of the preceding sections, we are now in a 
position to write down the Lagrangian density (22) in terms of the pseudo- 
complex valued heat amplitude, i.e., 

>C = ^(05*0* - 5*#*) - ^V0Vr. (43) 

Note that the Lagrangian density is real-valued. For later developments, it 
will prove useful to consider a slightly more general form, namely including a 
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quadratic term 

Lb = ^/(M0* - <fdt<^) - ^ V0V0* - 7#*, (44) 

which will be seen to give rise to an exponential clock with clock rate 7. Ob- 
serve that the Euler- Lagrange equations for (44) yield the equation of motion 
(9) for softly killed Brownian motion with a constant potential v{x) = 7. Then 
the canonically conjugate fields to 0, 0* are 



and we obtain the Hamiltonian density H. as the Legendre transform of £, 

n = TTdt(p + TT*dt(p* -C = Jv0V0* + 700*. (47) 

We observe that the Hamiltonian density is real-valued, such that the Hamil- 
tonian 

ii = \ d'^xH (48) 
can be consistently interpreted as the energy of the heat amplitude field. 



4 Single pairticle system 



For the intended development of the multi-particle theory in section 5, a thor- 
ough understanding of the one-particle system is inevitable. It is convenient 
to introduce a basis free Dirac notation for the heat amplitude equation (27), 
in terms of a Hilbert module formalism. 



^.1 Hilbert modules overP 

Let (7^, (■|-)]g) be a Hilbert vector space over IR. We define the associated 
Eilheit W-module {H, {■]■)) by 

n^TZa+eTZa-, (49) 

with an indefinite inner product (•!•) on 7i defined by 

{f\g)^{f_\g+)^a+ + {U\g.)^a^. (50) 



13 



It follows immediately that (-j-) is pseudo-hermitian and semi-linear in the 
first argument, 



{f\9)-{9\fr, (51) 
(A/|//5) = A>(/b), A,//eP. (52) 

As any real Hilbert space 71 has a basis, so has 7i, i.e., 7i is a free module. 
An operator A — A+a+ + A_a^ on Ji is seen to be self-adjoint with respect 
to (•!•) if, and only if, its zero-divisor components A± : TZ — > TZ satisfy 

= At, (53) 

where # denotes the adjoint with respect to the inner product on TZ. We 
further call an operator A : H — )• H real ii A — A*. It follows immediately 
that a self-adjoint real operator on 7i is a self-adjoint operator on 71. The 
spectral theorem for the Hilbert vector space 7Z therefore guarantees that the 
Hilbert module 7i decomposes into a direct sum of pseudo- complex eigenspaces 
Ex of any (i) sclfadjoint and (ii) real operator A on 7i, 

n^^E,{A). (54) 

Note that this slightly extended form of the spectral theorem for Hilbert mod- 
ules over P does not recur to the underlying real Hilbert space {7Z, (■|-)]{j) any 
more, but is entirely formulated in terms of the Hilbert module {7i,{-\-)). 
An important example of a Hilbert module is the space of square-integrable 
pseudo-complex valued functions on R*^, with the indefinite inner product 
{f\g) = / d'*x/*(x)g'(x). The operator X = x is self-adjoint and real, and 
hence gives rise to a complete basis of 7i. In contrast, the operator P = — /Vx 
is only self-adjoint, but not real. Hence the spectral theorem does not apply 
in the case of P. In the next section, we will see that the spectrum of P is 
indeed empty. 

Now let { |x)} be a complete orthonormal basis of eigenvectors of a self-adjoint 
real operator X : 7i — > 7i, i.e.. 



X|x)=x|x), (55) 
(x|y) = <5(x-y), (56) 

id-^^ j dx |x) (x| . (57) 

We will refer to components (x|/) of general states |/) G 7i with respect to 
the |x)-basis as the position space representation of |/). Stationary solutions 
(39) of the heat amplitude equation are elements |$) E H. Define states |p) 
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by their components in the position space representation, 

(x|p) = cos(px) + /sin(px), (58) 
and observe that we can, in turn, express the {|x)} in the |p) basis, 

|x) = 1 dp I - p) (p|x) . (59) 

Hence {|p)} constitutes another complete basis of H, which one may call 
twisted orthonormal, because of the relation 

(p|-q)=5(p-q). (60) 

It is sometimes useful to insert the identity in the form of the twisted or- 
thonormal basis, 

id^ = J dp \-p){p\. (61) 
4-2 Single particle Hilbert formalism 



We now put the abstract developments of the previous subsection to use in the 
discussion of a single Brownian particle. Define an operator P by its position 
space representation 

(x|P|y) = -/5(x-y)Vx, (62) 

and note its self-adjointness P"*" = P. The unitary operator exp(/a.P), acting 
on the Hilbert module, translates states in position space by a. We therefore 
identify its generator P as the momentum operator for the one-particle states. 
We observe that the states |p) satisfy the stationary heat amphtude equation 

H |p) = |p' |p) , (63) 
with H = ip^, such that the full time-dependent solution |p) is given by 

|p(t)) ^ eM~P't) |P) (64) 

Hence, any linear combination \(f){t)) = J dp(f){p) |p(t)) satisfies the time- 
dependent heat amplitude equation 

Idt m)) = H |0(t)) . (65) 

We emphasize that the states |p) are not eigenstates of the momentum oper- 
ator P, because we have 

P |p) = -P I - P) • (66) 
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This shows that the fundamental sohitions of the heat amphtude equation have 
no well-defined momentum, but a well-defined energy. Because the spectrum 
of — /Vx lies in i/M, there are also no linear combinations over with a well- 
defined momentum. This corresponds to the classical result that Brownian 
paths are almost surely nowhere differentiable. The probabihty density (3) for 
a Brownian particle to propagate from x to y within time t, is given in Dirac 
notation as the overlap of the initial and final states: 



after projection r+ to the heat part, according to our discussion of the heat 
projection operator r_|_ in section 3.2. Using that r+ : P — > IR is a ring homo- 
morphism, it is easy to show that the transition probabilities defined by (67) 
satisfy the Markov property (7) . This concludes our discussion of one-particle 
Brownian motion in terms of the heat amplitude equation. 



5 Multi-Particle systems 

We now want to give the idea of a multi-particle system a rigorous meaning. 
The basis for the understanding of a multi-particle system is the understanding 
of the one-particle system given by a Hilbert module (7i, (-I-)) and an equation 
of motion Idt \(p{t)) = H \4>{t)), as discussed in the previous section. In the 
one-particle theory, the state of a system is given as an element of the space 

n. 



5.1 Fock space 

The one-particle space Ti is now extended, in order to accommodate the case 
of zero, one, two, ... particles, the appropriate structure is the so-called Fock 
space generated by 7i 



where Ti.^ = P. We will consider indistinguishable particles obeying Bose 
statistics, i.e., an exchange of two particles goes unnoticed. This choice is 
reflected in considering the boson Fock space 




(67) 





n=0 i=l 



oc n 



(69) 



n=0 i=l 
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where denotes a symmetrized tensor product. The symmetrization removes 
any order of the factors in any of the spaces which therefore encodes 

the indistinguishabihty of the particles represented by these factors. Antisym- 
metrization would also lead to indistinguishable particles, but an exchange of 
an even number of them would not go completely unnoticed, but rather effect 
a sign change. 

The one-dimensional space HP accommodates a system of zero particles, and 
we call its normalized basis vector 

|0) = 1©O®O0O® ••• e 07i (70) 

the vacuum of the multi-particle theory. The vacuum is a new concept that 
does not occur in the one-particle theory. The boson Fock space 07Y inherits 
the inner product on 7i in a natural way, 

(Mi • • • Un\Vi Q---(DVn) = {ui\vi) . . . {Un\Vn) ■ (71) 

Define a set of linear operators acting on the Fock space 07i, labelled by 
elements of the one-particle Hilbert space Ti,. Their action on the n-particle 
subspace C Q7i is defined by 

4p> I""!) © • • • I""") = IP) kl) © • • • \Un) ■ (72) 

In particular, we can generate the entire one-particle space Ti from the vacuum, 

«tp) 1^^) = |P) • (73) 

The operators are called creation operators. Similarly, we define linear 
annihilation operators : QH — QH, labelled by elements of the dual 
space H*. They act on an element of the (n-|- l)-particle subspace 7Y®("+^) by 

n 

«<p| ko) • • • \Un) = Yl (Pl^i) l^o) © ■ ■ ■ \Ui) • • • \Un) , (74) 

where the hat on \ui) denotes the omission of this vector in the tensor product. 
In particular, all a" annihilate the vacuum, 

a^pl |Q) = ® ® ® . . . . (75) 

Linearity allows to extend the definitions of a+,a^ to the entire Fock space 
[15]. One easily checks that 

= (p|q) (76) 

from definitions (72) and (74). The choice of Bose statistics for the multi- 
particle spaces immediately requires 



."(p|'"|q>. 
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«{p| ' «(q| 



= 0, 
= 0, 



(77) 

(78) 



such that a^^a^^ |^^) = |p) |q) is well-defined. An important property 
for calculations is that creation and annihilation operators are adjoint with 
respect to the inner product (71), 



X 



for X e y G Ti®*^*^"*"^). Finally, we define an operator- valued field 



(x|p), 



(27r)'^/2 <Pl 

from which we can extract the annihilation operators by 

%| =(27r)'^/' /rfx $(x) (-p|x). 



(79) 



(80) 



(81) 



Note that up to this point, the creation and annihilation operator gymnastics 
have no physical meaning. They just describe the 'kinematics' of a multi- 
particle system, independent of any particular dynamics. 



5.2 Second quantization of heat 



We now want to derive the stationary dynamics of the free multi-particle 
theory from the stationary dynamics (63) of the one-particle sector H. To this 
end we require that the one-particle states (73) of the multi-particle theory 
still satisfy the equation (63), 

Hal) 1^) = IpSp) I") • (82) 

The same way in which we understand the onc-particle states in terms of 
creation and annihilation operators, we can construct the appropriate 'Hamil- 
tonian' H in terms of creation and annihilation operators, 

H^l c^'Py 4-p)«(p| , (83) 

such that (82) holds, as is easily checked using the commutation relations for 
the a+, a^. Note the sign of the label of the creation operator in (83), due to 
the pseudo-complex structure. It is now possible to express the Hamiltonian 
in terms of the operator- valued field 



H 



J d'^x^Vl'(x)Vl'+(x), (84) 
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where denotes the adjoint of $ with respect to the inner product on H. 

We know that the time-dependent one-particle states evolve according to (64). 
In the multi-particle theory, it is convenient to transfer this time dependence 
of the states to the operators a"*", a" or, equivalently, $. Clearly, this change 
to the so-called Heisenberg picture must preserve the inner product between 
states. If we denote the time evolution operator U{t) = exp{IHt), and note 
its unitarity UU^ = 1 with respect to the inner product on H, then the 
Heisenberg picture 0(t,x) of the stationary field $(x) must be given by 

0(x)=0(t,x) = t/(t)$(x)t/+(t), (85) 

in order to preserve the inner product on Ti. This definition can be rephrased 
as the Heisenberg equation of motion 

Idt^ix) = [h, 4>ix)] , (86) 

using that H is self-adjoint. Using the identity [A,BC] = [A, B]C + B[A,C] 
and (84), one can verify that the operator- valued field 4>{x) satisfies 

IdtU^) = ^A0(x). (87) 

This is the most important result of this section: We have shown that the dy- 
namics of the onc-particlc sector fully determine the dynamics of the operator- 
valued field 0(a;), i.e. the dynamics of the creation and annihilation operators! 
Indeed, the Lagrangian density for (87), 

C = ^(0+9*0 - 09*0+) - |V0+V0, (88) 

gives rise to the Hamiltonian (84). Note that now the ordering of the fields 
in (88) is relevant, and must be chosen as above to yield (84). An exphcit 
integral representation for the Heisenberg picture field (f){x), whose equations 
of motion follow from the more general Lagrangian density (44), including a 
quadratic term 700+, can be obtained by observing that 

^X = «p(^-^P'-7r, (89) 

which is immediate by an induction argument. It follows that for these slightly 
more general dynamics, the Heisenberg picture field is given by the explicit 
expression 

H^) = J ^^^j^apexp(-/(-p^ + 7)a;°) {cos(px) + /sin(px)} . (90) 

For completeness, we remark that the field and its canonically conjugate 
field TT = -Sfr = — ^0+ satisfy the equai time commutation relations 

dat4> ^ 



19 



[0(t,x),7r(t,y)] = -^(x|y), (91) 

[0(t,x),0(t,y)] = O, (92) 
[7r(t,x),7r(t,y)] = 0, (93) 

as is easily checked from the commutation relations of the time-independent 
a"*" and a~ operators. These equal time commutation relations contain no 
dynamical information, but simply encode the multi-particle kinematics. The 
next subsection deals with the explicit calculation of the commutator of fields 
at different times, yielding dynamical information on the system. 



5.3 Propagator for free Brownian motion 



We are now in a position to determine the one-particle propagator for Brown- 
ian particles within the formalism of the multi-particle theory. The result will 
present a check on our construction, as it must coincide with the transition 
amplitude (67) obtained in the one-particle theory. However, the calculation 
demonstrates the abstract concepts developed in the previous section at work, 
and we therefore present it in some detail. In order to slim down the notation, 
we will omit the hat on the operator-valued field from now on. Wc are no 
longer dealing with the classical heat amplitude 0, so that there is no danger 
of confusion. Consider the amphtude for a particle to be created at spacetime 
point X and to be destroyed at y, 

B{x,y) = {n\(t>{y)(t>+{x) \n) . (94) 

It is convenient to rewrite this in terms of a commutator 

B{x,y) = {n\[(t>{x)A^{y)]m (95) 
for the Heisenberg field 0. Prom the expansion (90) we obtain 

B{x,y) = j ^exp(-/(lp2 + 7)Ai)cos(p.Ax), ^gg^ 

where At = y° — x^ and Ax = y — x. We are interested in the transition 
amplitude for a Brownian particle, and hence must project B{x,y) by F^, 
yielding 

r+B{x,y) = J ^exp(-(ip2 + 7)At)exp(zp.Ax), (97) 

where we have used the symmetry of the integral in p to extend the cos (p. Ax) 
to exp(p.Ax). Completing the square and shifting of the integration variable 
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P^P + i^ gives 



(Ax)2 



2At 



/ 




(98) 



Rd. {2'kY 



We define the retarded propagator 



BR{x,y) = 9{y' 



,0 



x')B{x,y), 



(99) 



with ^(0) = 0, and extract the transition probabihty from the retarded prop- 
agator via the projection (37), 



This is exactly the transition density for Brownian motion in d spatial dimen- 
sions for 7 = 0, already obtained in the one-particle formalism. For 7 > 0, 
the Brownian particle is seen to carry an independent exponential clock with 
clock rate 7. Later, we will put this observation to use in the construction of 
branching Brownian motion, which possesses the Markov property if, and only 
if, the lifetime of each generic path is exponentially distributed (see section 



Also observe that the multiplicativity of r+ corresponds to the semigroup 
structure of Brownian transition operators. The Feynman propagator is de- 
fined as the time-ordered product of field operators. 



p{x,y) = T+{BR{x,y)) 




(100) 



7). 




(101) 



Using the integral representation 




(102) 



one finds 




(103) 



where k = (A;°, k) and the inner product is euclidean. 
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6 Interacting Multi-Particle Systems 



Our criticism in section 2.4, of the standard treatment of Brownian processes, 
was focused on the fact that such models are not fully specified by their 
dynamics, but also require an undesirable case-to-case adaptation of the un- 
derlying mathematical structures. The pseudo-complex Fock space formalism, 
developed in section 5, provides the solution to this problem. One can now 
study almost arbitrary extensions of the operator- valued field equation (87), 
in order to describe multi-particle systems, where the particle number changes 
dynamically. In the field theoretical framework, such dynamical creation and 
annihilation of particles are already fully encoded in so-called interaction terms 
added to the free Lagrangian. In this section, we first develop the theory a 
little for general interaction terms >Cint- example, we will study the bi- 

nary branching and dying model in the next section. 



6.1 Unitary time evolution 



Consider Brownian motion with a general interaction term ^i^t' giving rise to 
a total Lagrangian density 

£ = £B + Ant[0,0*]. (104) 

In terms of the Hamiltonian formalism, this gives rise to an additional inter- 
action term 

Hint = - j d^xCnt (105) 

in the complete Hamiltonian H = Hq + i^int- It is convenient to define the 
so-called interaction picture field 4>i as the Heisenberg picture field of the free 
theory, 

M^) = exp(7//o^)^(x) exp{-IHot), (106) 
where x = (T, x) Hence, the full Heisenberg field of the interacting theory can 
be expressed in terms of 0/ as 

(f)(x) = V{ty(j)i{x)V{t) (107) 

with 

V{t) = exp{IHot) exp{-IHt). (108) 
The unitary evolution operator V{t) satisfies the differential equation 

-IdtV{t)^Hi{t)Vit), (109) 

where 

Hj{t) = exp(7//oi)^int eM-IHot) (110) 
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is the interaction Hamiltonian in the interaction picture. Equation (109), with 
the initial condition \^(0) = 1 has the iterative solution 

V(t) = 1-1 f dtiHAtiWit). (Ill) 
By iteration one gets, with = t, 

OO ™ / ft —I \ 

V{t) = Y.i-irl[(l~ dUH:{U)), (112) 

n=0 i=l ^-^O ^ 

such that the introduction of a time ordering operator and appropriate com- 
binatorial factors [10] yields 



OO (_T)n 



/*° dt'Hi{t') 
Jo 

= rexp (^-I j'j dt'Hi{t')^ . 



(113) 



6.2 Correlation functions 



The field (f) encodes the full time-dependence of the interacting theory, whereas 
(pi simply describes the free theory, according to its definition (106). The 
(n, m)-point correlation function of the free theory is defined by 

(n| T{(j>i{x^) . . . MXn)(l>t{Xn+l) . . . 0|(^n+m)} |^^) , (114) 

such that the Feynman propagator (101) is recognized as the free two-point 
correlation function. The vacuum \Q) is the ground state of the free theory, i.e., 
it is the eigenstate of Hq with minimal eigenvalue. Correspondingly, we define 
the vacuum of the interacting theory, denoted by |E), as the lowest eigenvalue 
eigenstate of the full Hamiltonian H. Using (113), it is a standard exercise [10] 
to show that the (n, m)-point correlation function for the interacting theory, 

(S| T{0(a:i) . . . 0(0:00+ (xn+i) • • • 0+(^n+m)} |S) , (115) 

where is the Heisenbcrg picture field (107) of the interacting theory, can be 
calculated from the interaction picture field 0/ (106) and the free vacuum \ fl) 
as 

(Q| T {M^i) . . . 4>i{Xn)4>t{Xn+l) ■ ■ ■ 4>'t{Xn+m) exp(-/ / d''^^ zHljz))} \9) 

(Q| Teycp{-I ! d<^+^zni{z)) \Q) ' 

(116) 

where Tii is the Hamiltonian density for the interaction Hamiltonian Hj. It is 
desirable to simplify this expression by replacing the calculationally awkward 
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time ordering operator T by so-called normal-ordered expressions. The follow- 
ing subsection introduces Feynman diagrams as the appropriate technology. 



6.3 Wick's Theorem and Feynman diagrams 



We introduce the normal ordering operator A?", which acts on a product of 
creation and annihilation operators by regrouping the factors such that all 
annihilation operators are sent to the right of all creation operators, e.g, 

Wick's theorem states how a time-ordered product of heat amplitude fields 
0/, as occurring in (116), can be expressed by normal ordered expressions and 
contractions. 



T [(t)l{xi)(t)l{x2)(t)j{x^) . . . 



(f>l{xi)(f)i{x2)4>t{x3) . . . (l)l{Xf ) 



+N [all (possibly incomplete) contractions] , 
where a contraction is defined as 

C[Mx)cPt{y)]=BFix-y). 



(118) 
(119) 



(120) 



Note that for Brownian processes, only future directed contractions contribute 
in the above expression, because we find from (99) and (101) that Bp = B^. 
As an example, assume the following string of fields is already time-ordered, 
then 



T [<i>i{a)<i>j{h)<i>i{c)<i>j{d)} = N [0|(6)0+(rf)0,(a)0/(c 

+5f(a,6)Ar[0,(c)0+(rf)" 

+BF{a,d)N[(l>i{h)^(t>i{c) 

+BF(c,d)N[Ma)(l>t(b) 
+BF(a,b)BF{c, d). 



(121) 



A proof of Wick's theorem can be found in [10]. Representing the propagator 
Bf{x, y) by a directed line joining the points x and y, 

« = BF(x,y) 



we can write correlation functions in diagrammatic form. Note that the normal 
ordered terms will not contribute to the correlation function, as the sandwich- 
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ing between vacuum states annihilates them. The integral over Tij in (116) 
will lead to diagrams with vertices, where all occurring vertices are integrated 
over. 

Before we undertake the enterprise of calculating correlation functions for a 

particular theory, however, we can use Wick's theorem to simplify expression 
(116) significantly. As only future- directed contractions survive, we immedi- 
ately have 

{n\ Texp{-I J d'^+^zHiiz)) \Q) = 1, (122) 

and hence formula (116) for the computation of correlation functions of the 
interacting theory simplifies to 

(1^1 t|0,(xi)...0,(x„)cxp(-/ J d^'+'zHiiz))^ (123) 

This is consistent with a general result for any interacting quantum field the- 
ory, namely that the denominator in (116) cancels all disconnected vacuum 
diagrams in the numerator (see e.g. [10]). In the present case of Brownian 
motion, however, we saw above that the Feynman propagator is retarded, and 
so there simply are no non- vanishing vacuum bubbles. 

For general interactions Tii, it is non-trivial to calculate the (n, m)-point cor- 
relation functions, and very often it is impossible to do so exactly. In those 
cases, which present the rule rather than the exception in realistic quantum 
field theories, one must resort to a truncated series expansion of (123). The 
diagrammatic representation of a summand in such a perturbation series is 
called a Feynman diagram. 

In the next section, we will show that the binary branching and dying model 
is exactly solvable. Nevertheless, we will first expand the corresponding cor- 
relation functions, and then exactly sum the resulting series by integration of 
a set of recurrence relations. 



7 Example: Binary Branching Brownian Motion 

7.1 Feynman Rules 

Binary branching and dying processes of Brownian particles are easily for- 
mulated as interaction terms in the field theoretical treatment of Brownian 
motion. This model is completely specified by the Lagrangian density 

C = jCb[(I>] + + 7/30>0, (124) 
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where 7 G M.^, a,j3 E [0,1], and a + (3 = 1. The fertihty distribution 
{pk '■ k = 0,1,2}, encoding the probabihty that a dying particle will give 
birth to k children, will turn out to be p = {a, 0,(3). The parameter 7 will 
be identified below as the clock rate of the exponential clock controUing the 
process. Note, that the equation of motion corresponding to C is the reaction 
diffusion equation (14) with p = {a, 0, f3) for classical binary branching Brow- 
nian motion. The free theory Cb is as defined in (44). Applying the analysis 
of the previous section to the model (124), one obtains the position space 
Feynman rules [10] 

z 

as the fundamental building blocks of the trees generated by the branching 
process. Thus, at last, we recover the trees of branching Brownian motion 
from the field theoretical formulation! It is an interesting question of whether 
and how one can extract information on the genealogy of the process from the 
Feynman series, or indeed the exact summability shown below. The develop- 
ment of such ideas, however, is beyond the scope of the present paper. Before 
studying more complicated questions within the binary branching model, we 
calculate the probability that a single particle will either branch or die within 
a specified time interval. 

7.2 Exponential Clock 

The amplitude for a particle starting at spacetime point x to die within a time 
interval At is given by 

/7a dy^ j d'^yBpix, y) = Ia{l- e"^^^) , (125) 

using the Feynman rules above. Similarly, one obtains the probability for a 
branching event within the time interval At. Hence the total probability for 
the occurrence of branching or dying within At is found by the r+ projection 
as 

1 - e"^^^. (126) 
This exponential clock behaviour is crucial for the process being of Markov 
type, as explained in section 2. In the field theoretical treatment, it need not 
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be imposed, but is enforced, unless 7 = in (124), but then all interaction is 
switched off and one is left with free Brownian motion. Hence, the description 
of binary branching by the Lagrangian density (124), within the field theoret- 
ical framework, implicitly contains the branching tree, the exponential clock, 
and the fertility distribution, which must all be separately specified in the 
conventional approach. 



7.3 Exact summation of the binary branching Brownian motion 



The fully interacting model (124) provides a non-trivial test bed for our for- 
malism. Here, the one-point function already has contributions from any order 
in a and (3, 




where the subscript t denotes the earliest time at which all particles occurring 
in the series have died. Similarly, the two-point correlation function is 



+ + — + + 



with contributions of order C(a) = C(/?). 

We will put to use the special properties of the model at hand to exactly sum 
these Feyman series. The one-point function A{z) and the two-point function 
D{x,y) of the interacting theory satisfy the coupled recurrence relations 




Algebraically, these read 
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D{x,y)^F{x,y) + I-fP J^^ dw^ J d'^w BF{x,w)Ay^{w)D{w,y). 



The translational invariance of the theory imphes that the one-point function 
At{z) only depends on the difference r = i — 2;°, and hence we denote At{z) 
as A{t — 2;°). The recurrence relation for At can then be written 

A{t) = I'ja r dw^e-^""' + /7/3 /" (iw°e-^^°^2(T - (127) 

JQ Jo 

using that the spatial integral over Bf is normalized to unity. Multiplication of 
(127) by e'*"^ renders the second integrand independent of r. By differentiation 
with respect to r we obtain the differential form of the recurrence relation 

(127) . Writing A for the projection T^A, it reads 

= 7a-7i(r) + 7/?i2(^)^ i(0) = 0. (128) 

ar 

Indeed, equation (128) is equal to (16) for ^ = 0, because there we have 
W{ip) = a + (3(f'^ for a binary branching mechanism. The unique solution of 

(128) subject to the constraint a + P — 1 is 



1-:^' iff « = = ! 



(129) 

^ (1 - VI - 4a/:/tanh ^^/T^^Ao^'^t + C 



A{t) 



where tanh(C) = (1 — Aa(3)~^. 

From this result we can easily extract the extinction probability of the process, 
i.e., the probability for the process to have died completely at temporal infinity 
r — >• 00: 

_ l-|l-2a| I -r^ for a < i 

with almost sure extinction for a ^ ^ ^ as one would intuitively expect. 
Indeed, (130) is a classical result for branching Brownian motion [16]. As the 
function D{x, y) only depends on the difference y — xoi the space-time points 
X and y, we may write the recurrence relation for D{x, y) in the form 

D{z) = B{z) + 77/3 dw° J d'^wB{w)A{z^ - w^)D{z - w). (131) 

Requiring D to be integrable, standard arguments using Banach's fixpoint the- 
orem, guarantee existence and uniqueness of a solution to the integral equation 
(131). 
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7.^ Extension to Brownian catalysts 

The extension of the model to Brownian catalysts is straightforward in the field 
theoretical treatment. One simply provides terms for both the free Brownian 

particle and the free Brownian catalyst ip, say. The type of interaction 
between them is then modeled in the interaction terms. As an example, self- 
intoxicating particles are modeled by 

C = + Cb[iI^] + -far(t>* + l(^(t>*H. (132) 

where ip presents the toxic substance. The one-point function for this model 
is exactly solvable, now due to its simple two-point function, which can be 
exactly summed and immediately gives the one-point function. The analysis 
of catalytic models will, in general, be more involved than the exactly solv- 
able binary branching or self-intoxicating model. However, expressing catalytic 
Brownian motion in the form of a pseudo-complex field theory, most of the 
apparatus of quantum field theories is applicable to these questions. This is 
expected to allow insights into such processes that are hidden in the standard 
formalism. 



8 Conclusion 

The product structure of heat phase space emerges as the key to a geometrical 
understanding of characteristic properties of Brownian motion, most notably 
the non-differentiability of Brownian paths. This structure is concisely en- 
coded in the commutative ring of pseudo-complex numbers, giving rise to 
dynamics for a pseudo-complex valued heat amplitude. The real-valued heat 
can be extracted from the heat amplitude by a unique additive and multi- 
plicative projection, reflecting the absence of interference and Markov prop- 
erty for heat propagation. The discussion of an abstract Hilbert module over 
the pseudo-complex ring provides a convenient formulation of one-particle 
Brownian motion in terms of heat amplitudes. The appropriate kinematical 
framework for the discussion of Brownian processes with non-constant parti- 
cle number is the pseudo-complex Fock space generated by the one-particle 
system. The operator-valued second quantized field equations are rigorously 
constructed and extended by interactions. A standard derivation of the corre- 
lation functions for interacting theories in terms of the free fields yields their 
diagrammatical representations. 

The central result of the paper is the realization that multi-particle Brownian 
processes can be conveniently studied as pseudo-complex quantum field the- 
ories, with Brownian particles emerging as the quanta of the heat amplitude 
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field. Arbitrarily complex models are easily written down in the formalism. 
Indeed, the mere specification of the model dynamics, by a Lagrangian, gen- 
erates all elements of the conventional constructions: The Feynman diagrams 
of the pseudo-complex quantum field theory naturally generate the trees of the 
standard approach; interactions cannot be switched on without also starting 
an exponential clock, which guarantees the Markov property of the process; 
the fertility and death rates are concisely encoded in the dynamics of the 
model. No adaptations of the formalism must be made if the dynamics are 
extended by additional fields and arbitrary local interactions. This is in con- 
trast to the conventional approach discussed in section 2, and indeed solves 
the discussed shortcomings of the latter. 

The well-known binary branching Brownian motion provides an ideal test bed 
for the field theoretical description, and its one-point function is explicitly cal- 
culated, illustrating the extraction of information in the presented formalism. 
The extinction probability of binary branching Brownian motion is a point in 
case and yields, of course, the classical result. 

A field theoretical treatment easily allows the discussion of dynamical cata- 
lysts, which attract much interest in the present probabilistic literature on the 
subject [8]. The application of the theory to quantum systems immersed in a 
heat bath, a question that currently achieves a lot of attention in the context 
of quantum information theory, is under investigation. 

An important restriction of the systems that can be rigorously discussed, in the 
context of mathematical questions, is the requirement of exact solvability of 
the Dyson-Schwinger equations. However, for phenomenological applications 
of branching processes to practical problems, a diagrammatic expansion may 
just be the appropriate tool for extracting information on a complex system. 

Exerting some care in the application of standard field theoretical results to 
the pseudo-complex case, much of the vast apparatus of modern quantum field 
theory is now applicable to problems on Brownian motion, and vice versa. 
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